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^ ■ Abstract 

The depletion of gluons as the parton flux traverses a nucleus in a heavy-ion collision 
can influence the production rate of heavy-quark states. Thus the suppression of J /ijj 
lO ■ can be due to gluon depletion in the initial state in addition to nuclear and hadronic 

^^ \ absorption in the flnal state. A formalism is developed to describe the depletion effect. 

f~^ ' It is shown that, without constraints from other experimental facts beside the J/ijj 

C^ ■ suppression data in pA and AB collisions, it is not possible to determine the relative 

importance of depletion vs absorption. Possible relevance to the enhanced suppression 
seen in the Ph-Ph data is mentioned but not studied. 
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§ ; 1 Introduction 

The subject of J/ip suppression in heavy-ion collisions has been extensively investigated ever 
rS \ since its first proposal as a signature of color deconfinement [Q. The recent measurement 

c^ ■ of enhanced suppression in Ph-Ph collisions by NA50 has added considerable excitement 

to the possible interpretation of the data as such a signature |P, ^. Many alternative in- 
terpretations of the data have also been proposed |^, ||, |^ ||, 0. While some of them may 
have inconsistencies with all the pA and nuclear data |T^ , as pointed out in , a definitive 
interpretation of the Ph data has not yet reached general consensus. It is not the purpose of 
this paper to add to the controversy; in fact, the anomalous suppression in the Ph data is not 
our main concern here. We want to point out that there is a loophole in the interpretation of 
the pA and nuclear data (prior to the NA50 result) that is generally accepted, i. e. , the J /ip 
suppression is due to the absorption effects of the nuclear (and hadronic) matter that the 
cc system must pass through after it is produced. We investigate the possibility of another 
mechanism of J /if) suppression that has not been widely considered. It is the depletion of 
gluons before the formation of the cc state that leads to J /ip. If this new mechanism is found 
to be relevant to any heavy-ion collisions, including the NA38 experiments using O and 5* 



beams ITT], then the phenomenology of those past experiments must be re-examined before 



a definite conclusion can be reached concerning the anomalous suppression seen in the Pb 
data. 

The essential point to be made in this paper is that what happens to the gluons in the 
nuclei (apart from shadowing) before the basic subprocess g + g —>■ c + c is as important 
as what happens to the cc state after its formation. Most investigations on the subject 
concentrate on the latter, but the relevance of the former can easily be seen by considering the 
following extreme case. Suppose that an ordinary nucleus A collides with an extraordinary 
target nucleus B which is infinitely large. Clearly, the constituents of A cannot propagate 
through B indefinitely without momentum degradation and depletion. At some penetration 
depth the subprocess g + g —* c + c just cannot take place. Thus the production rate of cc 
(regardless of its fate afterwards) depends on the size of B and where the production points 
are. If that is accepted, then the issue becomes only a quantitative matter. What are the 
sizes of A and B when the initial-state effects are not negligible? 

There are two aspects of the initial-state effects on the partons: degradation and de- 
pletion. The degradation of parton momenta has been suggested previously [|12|, [1^. The 



mechanism of momentum loss relies on the radiation of soft gluons, as the partons pass by 
scattering centers. However, such processes of multiple emission of soft gluons take time and 



have been shown to be suppressed by the Landau-Pomeranchuk-Migdal effect ||Tj] , although 



the energy loss is not entirely negligible |T5[ . Indeed, the dependence of Drell-Yan production 
in p-A collision is oc A^ |l^ ; it may be taken as evidence of the ineffectiveness of multiple 
small-angle scatterings of quarks. Gluon depletion, on the other hand, is different. Whereas 
a quark undergoing scattering must remain as a quark, a gluon can, in addition to emitting 
gluons, also create qq pairs as it interacts with target partons. When that transmutation 
occurs, the gluon is lost from the beam, and the distribution of the gluons available for cc 
production downstream is thereby altered. For J/ip production the relevant momenta are 
> 1.5 GeV in the cc rest frame. If the subprocess oi g ^ qq involves an energy change of 
AE > 0.8 GeV, then the corresponding At (< 0.25 fm/c) is short enough for the subprocess 
to be completed in a distance corresponding to a mean free path A, i.e., in Az = A/7, where 
A ~ 2.7 fm, and 7 is the Lorentz factor (^ 10) for the CERN-SPS energy. Even ii AE is < 0.8 
GeV so that the formation time for qq is long, those gluons that produce the qq upstream 
cannot be effective in producing cc downstream in the same nucleus. For the dominant soft 
processes where the qq pairs are formed outside the nucleus, those quarks and antiquarks 
cannot contribute to the production of lepton pairs. Thus it is quite possible that the gluons 
can be stripped away from the incident gluon flux, as it traverses the target nucleus, without 
necessitating an enhancement of the dilepton production rate. Gluon depletion is therefore 
the loss of gluons from the incident beam along its path for any energy change, leading to a 
suppression of the g + g —>■ c + c subprocess downstream. 

We know that the conversion of gluons to sea quarks must take place efficiently in soft 
processes, since the gluons that carry roughly half the incident momentum in pp collisions are 
all turned into soft pions via the enhanced qq sea with the same total momentum fraction 
|T^, while the valence quarks produce the leading baryons, with no detectable glueballs 
produced. In pA collisions the wounded nucleon model that is successful in describing soft 
pion multiplicities |]rB| can be recast in the framework of the parton model, and one obtains 



a picture that is consistent with gluon depletion in that the incident gluons, once interacted, 
or converted to qq pairs, are ineffectual in producing more pions in subsequent collisions. 



The idea of gluon depletion was first applied to the problem of J /ip suppression in pA 
collisions by the use of an effective gluon distribution whose deviation from the gluon dis- 



tribution in the physical nucleon increases at larger A |T9[ . It is shown that the idea cannot 
be ruled out by the existing data on J /ip and T production rates in pA collisions. We now 
want to present a more detailed analysis of the problem for AB collisions. 

If gluon depletion is important in pA and AB collisions, then there should be a suppression 
of open charm production. There are some data on open-charm production, though not 
abundant. They are all on the single-inclusive production of the D mesons |20[-|23|, none 



on DD pair production. They reveal the nuclear dependences, A", that are characterized by 
values of a ranging from 0.81 to 1.02. The uncertainty is too large to be conclusive about 
gluon depletion. For comparison the observed A dependence for J /ip suppression corresponds 
only to a = 0.92 [Q. Furthermore, since the single-D inclusive production cross-section can 
include contribution from the hadronization of the cc component in the incident nucleons and 
from processes not initiated by gluon fusion, those data, even if accurate, cannot provide us 
with a reliable inference on gluon depletion. We therefore urge dedicated experimental efforts 
to examine the A dependence of two-particle back-to-back production of DD. Information 
acquired in such experiments can provide crucial constraints that can resolve some of the 
ambiguities uncovered in our study in this paper. 

2 Eikonalized Gluon Depletion 

The usual expression for the production of cross section of heavy-quark pairs QQ is 

^ = Y1 J dXidx2Fi{xi, Hf)Fj{x2, fiF)^ij{Xl,X2, flR) (1) 

where aij{xi,X2, ^r) is the cross section for the hard subprocess, i+j—^Q + Q,i and j 
being the partons involved and /ir being the renormalization scale. -Fj,j(xi 2,A*f) are the 
number densities of the partons i and j at momentum fractions Xi and X2 and factorization 
scale fip- Usually, the two scales are set equal to fiR = fip = fJ' = 2mQ. Equation (|lD, which 
is used for hadronic collisions has generally been applied to nuclear collisions also with the 
appropriate replacement of the parton distributions by Fi/A{xi) and Fj/b{x2) that take into 
account the shadowing effects in the nuclei A and B. Suppression of the detected onium 
states of QQ due to processes that take place after the production of QQ does not alter (|1|), 
which describes only the initial states of the hard subprocesses. 

The basic point about gluon depletion is to question the validity of (^ in the gluon sector. 
More specifically, the challenge is in the factorizability of the nuclear gluon distributions. If 
A and B are hypothetically large, then factorization cannot be valid on physical grounds. 
We give below a formulation of its nonfactorizability in terms of a physical cross section. 

Let us denote the nuclear thickness of nucleus A at impact parameter Ba by 



TA{bA)= dzpA{bA,z) , (2) 

where pA{bA,z) is the nuclear density, normalized such that 

jd%ATA{hA)=A . (3) 



For A > 1, we define 



TXibA, za)=[1-\] 1^ dzpAibA, z) , (4) 



A 



TX{bA,ZA) = [1- -r) dzpA{bA,z) , (5) 



A 



— oo 



so that / d'^bA [^a (^a, za) + T^(&a5 ^a)] = A — 1. The variable z is positive in the direction 
of A's momentum in the cm system. Thus if a hard subprocess occurs at za, then T^(6a, -^a) 
measures the nuclear matter in the path before the interaction point, while T^{bA, za) refers 
to the matter that trails behind. Clearly, the former is relevant to the initial-state interac- 
tion, and the latter the final-state interaction. Similar expressions are defined for Tg (bs, zb), 
except that zb is positive in the opposite direction, i. e. , in the direction of i?'s cm mo- 
mentum. Assuming that the hard subprocess is sufficiently rare so that in any AB collision 
it can occur at most once, we may identify za with zb at the interaction point, although 
the two variables are later independently integrated over to account for all possible relative 
positions in the A and B nuclei. 

The average number of inelastic collisions that a nucleon in A at (6^) ^a) suffers as it 
traverses i? in a straightline path leading to zb is ainT^{bB,ZB), where a^^ is the inelastic 
pN collision cross section. Strictly speaking, after the first collision of that nucleon (call it 
p) in A with a nucleon in B, the former becomes a broken nucleon (call it p') as it proceeds 
through the remaining part of B, and the relevant cross section for the subsequent collisions 
should be af^ rather than af^ ||2^. We do not make the distinction here now, and denote 
it by the generic symbol (Tin. Later, it will be combined with some other unknowns in the 
problem and become one overall parameter. 

The probability that p makes z/i collisions in B before getting to zb is 

71"z.i(&b,^b) = — r cri^Tj^{bB,ZB) exp -cri^Tj^{bB, zb) , (6) 

where a Poisson distribution has been assumed. After a collision the gluon distribution is 
modified by a factor h{xi), for which we assume a rather general form 

h{x,) = hoxUl - x,y . (7) 

Assuming that the same h{xi) applies at every collision, the overall modified gluon distribu- 
tion after z/i collisions is then 

G,A^i) = h''^{xi)gAixi) , (8) 

where ^^(^^i) is the gluon distribution of a nucleon in A before any collisions but with 
shadowing effects taken into account. Similar modification takes place for the gluons in B 
due to collisions with nucleons in A, and we have 

Gu2{x2) =h''^{x2)gB{x2) , (9) 



where V2 is the number of colhsions that a nucleon in B encounters in A before reaching za 
with probabihty ■n^.^ipA^ZA)- 

These modified gluon distributions are what must replace FiFj in (|l|), if gluon depletion 
is to be taken into account. Thus focusing on the gg —>■ QQ subprocess in (|I|) we have for a 
particular interaction point in A and B 



OO CXD 



aQQ{bA,ZA]bB,ZB) = J2 H ^'^i(^S'^-B)7r;.2(^A5^A) 

ui=0u2=0 

dxidx2G,,^{xi)G^^{x2)crgg^QQ{xi,X2) • (10) 

It is clear from this equation that factorization does not hold. The gluon density at xi (in 
A) depends on the path length in B from oo to zb, contained in TTi^^{bB,ZB)- The total 
production cross section of the QQ state is 



cr, 



QQ 



d%Ss dzA dzB Pa{s, za)pb{s- b, ZB)aQQ{s, za] s-b, zb) 



(11) 



where 6^ = s and bB = b — s. Equations (|TOD and (|TT]) represent an improvement of the 
initial-state description of AB collisions that has hitherto not been considered. 

For the absorptive effects on the production of J /if) we use the conventional description 
|p6| , ^ . First, let us replace o^gg^QQ{xi,X2) in (|^) by the cross section of the subprocess 
of J/ip production, agg^j/^{xi,X2), before any absorptive effect, but including intermediate 
states such as ccg, x, etc. |^. Next, we take the absorption into account by writing the 
survival probability in the exponential form 



exp|-cra TXibA.ZA) ^TB^bB.ZB) } 



(12) 



where Oa is the absorption cross section of J lt\} interacting with the final-state medium 
(whether hadronic, nuclear or quark-gluon plasma), leading to open charm. Putting all 
these factors together, we write the final result for J/^ production cross section in the 
following way: 



C"J/V 



d^bd^sdzAdzB Pa{s, za)pb{s — b, zb) 
dxi dx2 Fa{xi, s-b, zb)Fb{x2, s, ZA)cr, 



gg^J/'tp 



>i,a:2) 



■exp{-crii 



Ta{s,za) +Tj^{s-b,ZB] 



0"a 



T+{s,ZA)+T^is-b,ZB)]}.il3) 



where 



oo -I 

FA{xi,bB,ZB) = Xl ~r CrmTB{bB,ZB] 
1.1=0 ^1- 



i^i 



Gu^{xi) 



(14) 



oo -| 

FB{x2,bA,ZA) = J2 —i aiaTXibA,ZA] 
.2=0 '^s! L 



1^2 



Gu2{X2) 



(15) 



Equation (|T3), for example, can be given the interpretation of the (improperly normalized) 
gluon distribution of A modified by the depletion effects due to B. The extra normalization 



factor is the exponential term, which is now included in the last line of ([T3|) for a reason that 
will become self-evident in the next section. 



3 A Simplified Case at xf = 



In Eq. (|T3D we have an expression for the J/ip production cross section, obtained under a 
rather general description of the gluon depletion process in high-energy nuclear collisions. 
In principle, if all the parameters controlling different factors in the problem are known, the 
computation of aj/^ in accordance to (|13|) is straightforward. That is especially true if one 
determines only the suppression effects by computing the ratio 



qAB 



AB lABia) 



(16) 



AB(0) . 



where <J^u^"' is the J/ijj production cross section in AB collisions without absorption or deple- 
tion, since in that case the inaccuracies in the leading-order approximation of crgg^j/^{xi, X2) 
cancel in the ratio. 

At this point our investigation of the gluon depletion effects is still preliminary, since 
the various factors involved in the relevant dynamics are poorly understood. A systematic 
program for its exploration should therefore begin with a simplified calculation that can 
make transparent the connections between the physics issues and their phenomenological 
consequences. More detailed calculations can come later when proper focuses can be placed 
on specific issues, after a general picture becomes clear. Our immediate aim is therefore to 
capture that general picture and see whether gluon depletion can be relevant to the present 
and forthcoming experiments in the first place. 

The first step in our simplification is to consider J/ip in a narrow region around xp = 0. 
For y/s ~ 20 GeV that means xi ~ X2 — Mj/^/^/s ^ 0.15 or slightly higher for the 
production of cc state that can lead to J/ip by soft gluon emission. In the approximation 
that the integrations over xi and X2 in (jT^) need only be extended over the narrow range 
between 0.15 and, say, 0.18, beyond which d-gg^j/^{xi,X2) is negligible, we may replace the 
integrals by evaluating the integrand at Xi = X2 = 0.16, and obtain in consequence of (8) 
and (9) 



dxidx2G^^{xi)G,,^{x2)crgg^j/^{xi,X2) ^ cD 



U1+U2 



(17) 



where c and D are some constants. Actually c can depend on the nucleon numbers A and 



B on account of nuclear shadowing, but it will be canceled in the ratio Sjf^ . On the other 
hand, D represents the effect of gluon depletion and is raised to the power z/i + 1^2, thereby 
contributing a factor of crucial importance to us. 
Using (p!7D and ignoring c, we have 



/ dxi dx2 Fa{xi, hs, zb) Fb{x2, Ba, ZA)cr 



: exp jo-inL' 



'"gg^J/-(/'(^l?^2) 
<J,nTB{hB,ZB)\ [ai„T^(6A,^A)] /^"^+"^ 



T, 



^ ^ ... ^iA{bA,ZA)+TB{hB,ZB)\} . (18) 

Combining this simple result with the exponential factor in the integrand in (|13|) yields the 
probability factor 



P = 



= exp {-ad T4 {s, za) + Tj^{s- b, zb] 



-o-„ 



T 



''A{s,ZA) + T^{s-h, 



zb)]} 



(19) 



where 



(^d = 0-in(l - D) 



(20) 



We may call this the depletion cross section, inasmuch as 0"^ is called the absorption cross 
section. Note that a^ is an overall parameter, summarizing a number of imprecisely known 
factors. Among them is D. If there is no gluon depletion, then in (|^) ho would be 1 and 
a = P = 0. In that case we would have D = 1 and 0"^ = 0. In past investigations of J/ip 
suppression it is universally assumed that cr^ = 0. We now see how nonvanishing values of 
a and (3 can have phenomenological consequences. 

For phenomenology it is sufficient at first to start with (p!9D , which involves just two 
parameters, ad and aa- The thickness functions T^^ accompanying ad and aa are just what 
they should be. T^ is the nuclear thickness of A before the interaction point, and T^ is that 
for B, while T^b ^^^^ ^^le respective thicknesses after the interaction point. The symmetry 
of the suppression mechanisms is now complete: depletion during the pre-interaction phase 
and absorption during the post-interaction phase. Without further study it is not obvious 
which is more important. 

Integration over the geometrical variables can be significantly simplified without much 
sacrifice in accuracy, if we approximate the nuclear density by a constant value po, i.e., 
PA,B{r) = Po^{Ra,b — ^) where Ra and Rb are the radii of A and B, respectively. In that 
approximation the nuclear thicknesses are 



T^=p'^iLA±ZA) 

and similarly for Tg , where 



p^ = (i-iM)po 



(21) 



{Rj 



.2U/2 



^B 



{Rl- 



s-h 



vl/2 



(22) 



Combining ([13[) , (|T6D and (|l^) we have for the suppression factor 



SJf^ = NAhjd%dhU{h,s) 



(23) 



where 



[■La i-Lb r r n 

/ dzA / dzB exp -ad Pa {La - za) + Pb {Lb - zb 
J-La J-Lb ^ '- 

{La + Za) + pI {Lb + zb)] } 



p'a 



f^-2aap\LA _ ^-2aap\LA\ f^-2aaP%LB _ ^-2aaP%LB\ j 



P\P^B {(^d - (Taf 



(24) 



Nab = 4: f d%d'^sLa{s)LB{\s-b\) 



(25) 



The symmetry of the problem under the interchange of aa and ad is now explicit. Further- 
more, if ad = 0, (^3|) and (|2^) agree with the corresponding formula derived by Gerschel and 



Hiifner 261. 



Since we know very little about the dynamics of gluon depletion, we have no reliable 
information on the magnitude of ad- However, we do know that the Gerschel and Hiifner 
formula can fit the heavy-ion data on J /tp suppression, excluding the Ph-Ph collision result 
0, by use of Ua = 6-7 mb (and, of course ad = 0) [|]. We therefore can expect that when 
the effects of gluon depletion are considered, the combined suppression mechanisms would 
have roughly the same overall cross section. Nevertheless, we use the combined cross section 
defined by 

o-c = (Xa + ad (26) 

as a free parameter to fit the pre-Pb data. The ratio tj = ad/aa can still vary between 
and 1. The range r] > 1 leads to no new result because of the cTq <-^ ad symmetry of U{b, s). 
Thus the measurement of Sf,^ cannot resolve this ambiguity at this level of consideration. 
For definiteness we shall examine the range < r^ < 1. In that range it is not obvious by 



inspecting (p^)-([25|) how Sj,^ depends on A and B. A numerical computation is therefore 
necessary. 

4 Some Numerical Results 

A crude but quick estimate of (p3|) without doing the integrations is to replace La and Lb 
in the integrands by their averages, jRa and jRb, respectively pH], and p^s by Po for ^ 



and B > 1. Denoting the resultant approximation of SfPu by Sf,^, one obtains the analytic 
form 

where Xa and A^ are the mean free paths 

K = (o-aPo)"\ Ad = {adPoy^ . (28) 



For pA collisions (^) becomes 



2 



Using Ra = 1.2 A^/^ fm and p^^ = |7r(1.2)^ fm^, (P?]) and (P^) can be calculated as a function 
of A and B for various values of ?7 = ad/aa subject to the constraint (p6|). The numerical 
result for that will be given and discussed below. 

A direct computation of SfP^ by carrying out the integrations in ( pB]) is, of course, 
straightforward. Since the measurable total transverse energy Et depends on b, the integra- 
tion in b can be suspended by plotting the suppression factor against the mean longitudinal 
length L{b) at fixed b, which in turn can be related to the Ej-. L{b) is defined by 

Lib) = (La + za + Lb + zb),^,^^,^ 

/o^^ sdsJ^''deLA{s)LB{b- s) \La{s) + LB{b- s) 



J,''^sdsCdeLA{s)LB{b-s} 



(30) 



where 



:^(b-g) = (Rl-b'^- s'^ + 2bs cos eV^^QiRB- b-s\) , (31) 



G being the step function. Parenthetically, we remark that the value mentioned earlier for 
the average La is 

(L^)s,., = {I^Ra) J d's 1^^^ dzA [LAis) + za] = Ira . (32) 



Thus the approximation made in (pTl) is effectively 



im), ^ (La),,,,, + (Lb),,,,, • (33) 

The calculated results on Sf,^ obtained from the use of (27) or (29), as the case may 
be, and on the more precise Sfi^ from (23)-(25) turn out to be very nearly the same with 
differences being at the level of < 1%. We therefore present a single figure to represent both. 
It is shown in Fig. 1. The triangles are the calculated points for cXc = 6 and 7 mb, and 
1] = aa/cTa = 0.05. As expected, they compare very well with the data as shown in open 
squares, except for the Pb-Pb point. The solid straight lines are best fits of the triangles, 
and provide adequate fits of the data up to SU. That is the known result from earlier work 
before the Pb data, and is based on no gluon depletion. Now we fix ac at 7 mb and increase 
?7 to 1. The result is shown by the circles, which are fitted by the dotted line according to 
a quadratic formula of the form a + bx + cx^. Although the result does not differ too much 
from the f] = 0.05 case, there is nevertheless a perceptible bend downward at higher values 
of AB. We find that to be a very encouraging sign for the possible interpretation of the Pb 
data as being a manifestation of the gluon depletion effect, provided that an enhancement 
at high values of AB can be incorporated in an improved description of the effect. 

In Fig. 2 we show the results for Sf,^ calculated at specific values of b, but are plotted 
against L{b) using (30). The same values of cXc and r] as in Fig. 1 are used. The general 
agreement with the data is again very similar to that in Fig. 1, as it should. The fact that 
the f] = 1 points get within the error bars of two of the three Pb-Pb data points provides 
further motivation to take the possibility of gluon depletion seriously. 

The most striking feature in those Figs. 1 and 2 is that the dependences on r] are very 
small. Without further physics inputs it is not possible to extract from the data the realistic 
value of rj. A conservative conclusion is therefore that the possible contribution of gluon 
depletion to J ftp suppression cannot be ruled out. This statment is already of considerable 
importance in our view because firstly it means that a loophole has been found in the 
conventional approach to J/ip suppression and secondly with a crack opened in this new way 
of considering the suppression mechanisms it is possible to imagine enhanced suppression 
at large A and B, as will be discussed in the following section. The downturn of the Pb 
data requires an increase of either aa or ad- Any speculation on the increase of ad is not the 
purpose of this paper. However, until that possibility is firmly ruled out, the increase of o"a 
as the explanation of the Pb data should only be held as tentative, albeit a very attractive 
one. 



5 Comments 

We have raised the issue of gluon depletion in heavy-ion colhsions, developed a formalism 
to describe its effects on J/ip production, performed numerical computation to examine its 
consequences, and shown that the present data cannot exclude its possible contribution to 
J/ip suppression. The combined cross section is found to be ctc — 7 mb, but the ratio 
Tj = <7d/(Ja is undetermined because of the insensitivity of the suppression factor to t]. 

An independent experimental constraint is necessary to determine rj. We suggest the A 
dependence of back-to-back correlated production of DD in pA collisions. For DD produc- 
tion near threshold the formalism that we have developed is applicable, if (Jq is set to zero. 
Thus any observation of ctdi) oc A" with a < 1 would be a signature of gluon depletion. 

It is also possible to get extra information from J/ip suppression if we examine the xp 
dependence. The restriciton to xf = in Sec. 3 simplifies the problem so that the suppression 
factor Sf,^ does not depend on the detailed gluon distributions. However, for xp ^ 0, the 
details of all the factors in (p!oD will become relevant, and the effect of gluon depletion cannot 
be described by one collective parameter a^- Thus the subject has the potential of developing 
into a fertile field of phenomenology. 

A more pressing question is perhaps inescapable: does gluon depletion have any relevance 
to the more urgent issue of enhanced suppression observed in NA50 0? We have avoided 
addressing that issue in order to be clear about what constraints the pre-NA50 data can place 
on gluon depletion, the pertinence of which in heavy-ion collisions should be investigated 
independent of the NA50 data. We now ask whether there is any chance that the downturn 
of the suppression factor in the Pb-Pb data can be due to some aspects of the gluon depletion 
process. The formalism described in Sees. 2 and 3 does not lead to any prominent nonlinear 
behaviors in Figs. 1 and 2. However, it should be noted that the modified gluon distribution 
in (H) is obtained under the assumption that the same depletion factor h{xi) applies at each 
of the Ui collisions. That is a reasonable first-try to estimate the effect of multiple collisions, 
but it does not follow from any careful dynamical consideration. In the absence of a workable 
nonperturbative QCD calculation one can envisage a study in which gluon depletion is viewed 
as a gain-loss evolution process where gluons in a cell of momentum fraction xi ~ 0.15 are 
lost from the cell due to the g —^ qq subprocess, but the gain comes from higher-xi cells 
due to g ^ gg, for example. Since the initial gluon distribution gA^xi) behaves roughly as 
(1 — Xi)^, there are less gains than losses as ui increases. Thus it is conceivable that the 
modification factor in front of gAi^i) in (^ may increase as a nonlinear power of ui, resulting 
in an enhanced effect at higher Ui. The overall suppression when combined with absorption 
may show a break from linearity if below the crossover point the absorption is dominant, 
while above it the depletion is more important. These are speculations that need concrete 
calculations to gain substance. Nevertheless, unless such possibilities are excluded, there is 
no water-tight argument in favor of any other interpretation of the enhanced suppression. 

Despite the alluring challenges posed by the NA50 data, we feel that at this point it is 
more important to pin down the extent of gluon depletion in a systematic way, with emphases 
on pA collisions, correlated DD production, and the longitudinal momentum dependence of 
J/ip production. 
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Figure Captions 

Fig. 1 The suppression factor Sf,^, abbreviated as 5*, is plotted against AB for various 
combinations of the values of the combined cross-section ac and the ratio rj = cxd/c^a- 
The data are from ^. 

Fig. 2 Same as in Fig. 1, but plotted against L{b). 
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